On the singular spectrum of the Almost Mathieu operator. 
Arithmetics and Cantor spectra of integrable models. 



On ■ 

on • 
on : 

m ■ 
in ■ 



(N 

> 
m 
oo 
O 

ON 

o 

On 
On 



Oh! 

CD . 



X 



P. B. Wiegmann 

James Franck Institute and Enrico Fermi Institute of the University of Chicago, 5640 South Ellis Avenue, 

Chicago, Illinois 60637 
and Landau Institute for Theoretical Physics 

I review a recent progress towards solution of the Almost Mathieu equation (A.G. Abanov, J.C. 
Talstra, P.B. Wiegmann, Nucl. Phys. B 525, 571 ,1998), known also as Harper's equation or Azbel- 
Hofstadter problem. The spectrum of this equation is known to be a pure singular continuum with 
a rich hierarchical structure. Few years ago it has been found that the almost Mathieu operator is 
integrable. An asymptotic solution of this operator became possible due analysis the Bethe Ansatz 
equations. 



PACS number(s): 05.45. +b, 71.23.Ft, 71.30. +h 

1. Introduction In this lecture I review a recent 
progress |l[ towards decoding of one of the most puzzling 
strange set generated by a quasiperiodic Schrodinger op- 
erator : 



ipn+i + 4>n-i + 2Acos(# + 2TTni])ip n — Eip ri 



(1) 



The history of this equation as well as its applications in 
different branches of physics and mathematics are rich. 
This equation, known as Harper's equation, describes 
the electronic spectrum of one-dimensional quasicrystal 
(a particle in a quasiperiodic potential) and often used 
to study localization-delocalization transition (see e.g. 
P,|l). It also describes a Bloch particle in a uniform 
magnetic field and also known as Azbel-Hofstadter prob- 
lem ||,[|. It is a standard example of the operator, also 
known as almost Mathieu operator with a singular con- 
tinuum spectrum |^,^| . This list may be continued. 

The spectrum of this equation is complex. In the com- 
mensurate case, when rj is a rational number 

v = P/Q, 

one may impose the Bloch condition: 

if>„ = e lk ip n+Q 

Then, the spectrum consists of Q bands, separated by 
Q — 1 gaps. In the incommensurate limit, when 77 is an 
irrational number^] (P — > 00, Q — > 00) the spectrum be- 
comes an infinite Cantor se|j] with total bandwidth 
(Lebesgue measure of the spectrum) 4|A — 1 f|,|lll. 

The most interesting "critical" case appears at |A| = 1. 
Then the spectrum becomes a purely singular continuum 



'Although some properties of the spectrum depend on the 
type of irrational number n, here we consider typically Dio- 
phantine numbers. These numbers have a full Lebesgue mea- 
sure and thus sufficient for almost all physical applications. 

Closed, nowhere dense set with no isolated points 



0. In this case wave functions lost their extended char- 
acter and not yet localized but exhibit a power law scal- 
ing. Moreover, there is numerical evidence and almost 
a consensus, that in this case (A = 1) the spectrum and 
wave functions are mult ifr act al flPfl . 

Multifractal sets exhibit a sort of conformal invariancc 
and are expected to be described by methods of confor- 
mal field theory. This theory, however, is yet to be devel- 
oped and scaling properties of sets generated by dynam- 
ical systems and by closely related quasiperiodic systems 
are far from being understood. 

Since the empirical observations of Hofstadter || , the 
evidence has been mounting that the spectrum of (1) 
(Hofstadter butterfly) as well as other quasiperiodic equa- 
tions with a differential potentials are regular and univer- 
sal rather than erratic or "chaotic" . Few years ago it has 
been shown p|-|l5| that despite the complexity of the 
spectrum, the Harper-Azbel-Hofstadter-almost Mathieu 
operator ([!]) equation at any rational r\ = P/Q is in- 
tegrable and can be "solved" by employing methods of 
integrable systems Ansatz (BA). This had opened the 
possibility of describing the complex behavior of an in- 
commensurate system as a limit of a sequence of inte- 
grable models. I hope that this solution will help to apply 
conformal field theory to dynamical systems. 

The symmetries of the problem which eventually lead 
to its integrability are the most transparent it its "mag- 
netic" interpretation. Consider a particle on a two di- 
mensional square lattice in a magnetic field with a flux 
$ = 2ttt] per plaquette. Its Hamiltonian is: 

H = T x + T- 1 + \(T y + T~ 1 ), 

where operators T x and T y describe translations of the 
particle in x and y direction by a lattice site. In magnetic 
field translations form a Weyl pair: 



(2) 



where 



1 



The Harper's equation, then appears as a result of rep- 
resentation of translation operators as a shift and multi- 
plication: 



T x i> n = tpn+l 
Tylpn 



(3) 



2. Hierarchical tree. I begin by describing the scal- 
ing hypothesis for the spectrum of an incommensurate 
(quasiperiodic) operator with a purely singular contin- 
uum spectrum. To do so, we need a notion of the hierar- 
chical tree. 

Let us consider a sequence of rational approximants 
fj(i) — Pj/Qj with increasing Qj to an irrational flux rj 
so that: \rfi> — T]\ < c(Qj)~ 2 , where c is a j-independent 
constant^. A Harper equation taken for each iyi> is gen- 
erations of the hierarchy. Let us consider a graph (with 
no loop), which connects the fc-th band of the genera- 
tion rj 3 (the daughter generation) to a certain band k' of 
some previous (parent) generation r)^~ l K We call it a hi- 
erarchical tree if energies Ej (J7) belonging to any branch 
J of the tree form a sequence converging to the point 
E(J') of the spectrum in such a way that the sequence 
Q 2 ~ tJ \Ej{J) — E(J) \ is bounded but does not converge 
to zero. 

The set of numbers e j are anomalous exponents. In a 
multifractal spectrum, anomalous dimensions depend on 
the branch. They and the tree characterize ultrametric 
properties of the spectrum. 

Let us stress that the very existence of the hierarchical 
tree is a hypothesis and the tree constructed below is the 
conjecture. We call it scaling hypothesis. 

To construct the hierarchical tree it is necessary to find 
the sequence of generations and a rule to determine the 
parent generation and a parent band out of a given band 
of a given generation. In other words, the hierarchical 
tree is determined by a sequence of rational approximants 



■q and by a mapping 



(k,P,Q)^(k',P',Q'). 



(4) 



where k and k' are labels of a daughter's band and the 
parent's band of generations P/Q and P'/Q respectively. 
To describe the hierarchical tree we will need a notion of 
a discrete spectral flow and its rate Hall conductance. 

A heuristic definition of the spectral flow is as follows. 
Let us consider a spectrum of the problem with a given rj 
and choose some (big) gap. We label it by k. Now let us 
change rj by a small Srj, such that newly appeared gaps 
in the vicinity of the edge of the big gap, will be smaller 
than the big gap k. Then we can look on new levels 



"'"This sequence always exists and can be constructed e.g. 
from the Farey Series [jl6[ 



appeared within the big gap, close to the bottom of the 
gap. The number of these levels, i.e., the number of levels 
6Nk, crossing an energy E, lying inside the gap, close to 
its lower edge, is a spectral flow. One can also treat the 
spectral flow as a number of levels appeared within a 
"big" band adjacent to the gap from below. The rate of 
the spectral flow is 



5N k 
Si] 



(5) 



is known to be the Hall conductance of the gap p7| ] 
(Streda's formula) of two dimensional particles in mag- 
netic field. This number is an integer and depends on 
the gap. The difference between the Hall conductances 
of nearest gaps, i.e., the spectral flow into the k-th band 



a(k) = a k - a k +i 



(6) 



is the Hall conductance of the band. 

The index theorem identifies the Hall conductance of 
the band with the number of zeros of the wave function 
4> n (k, k') within the Brillouin zone: < k < 2ir, < k < 
2tt/Q, i.e., with the Chern class of a band: 



a(k) 



1 
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V In ?/>„(£;, k')dk, 



where the integral goes over the boundary of the mag- 
netic Brillouin zone of the k th band. 

The Hall conductivity of the k-th gap varies the range 
—Q/2 < o~k < Q/2 and obeys the Diophantine equation 
Plpl 



Pat = k (mod Q). 



(7) 



In it turns the Hall conductance of a band cr(fc) is allowed 
to have only two values. They can be found explicitly. 
Let us consider a continued fraction expansion of 



U) _ 



1 



[rii, n 2 ,n 3 , . .., Tij\ 



Then the Hall conductance of the gap k is 
Qj (■},{■ 1^ 



Cfc 



Qj 



-k ■ 1 



2 I ' 



(8) 



(9) 



while the two values of the Hall conductance of the band 
k may be 

o-(k) = {-iy- l Q^i, or (-lyiQi - Qj-i) (10) 

here {x} is fractional part of x. 

Now let us turn to the hierarchical tree. We conjecture 
that the hierarchical tree is the spectral hierarchy - an in- 
tegral version of the spectral flow. Let us consider two 
close rational P/Q and P'/Q' with Q' < Q. The number 



2 



of states per lattice site in a band is 1 /Q and l/Q' re- 
spectively. If there is a band k of the problem with a flux 
P/Q, such that its Hall conductance is a ratio between 
the difference of the number of states and the fluxes 



{Qi 



1. 



1,2,...} (15) 



11 P P' 



(11) 



then we say that the band k of the generation Pj Q has 
a "parent" band in the generation P'/Q'. The absolute 
value of the Hall conductance is the difference between 
number of states in the parent and a daughter band: 



Q-Q'= \a{k)\ 



(12) 



This formula may be viewed as an integrated Streda 
formula (|[). It determines the flux P'/Q' and and by 
virtue of an iterative procedure, generates a sequence 
of rational approximants (generations of the hierarchical 
tree), rfi> to r\. 

The integrated Streda formula is not enough to deter- 
mine the tree. We complete it by the adiabatic prin- 
ciple, which states that the levels do not cross each 
other along a tree. This proposition may be put in 
symbols. Let us enumerate all states from the bottom 
of the spectrum and characterize them by a fraction 
v = (number of the state)/Q. All states in the k-th band 
have the fraction (k — 1)/Q < v < k/Q. The adiabatic 
principal assumes that a state in the middle of the band 
k with a fraction (k — 1/2) /Q has a parent somewhere in 
the parent band k' , i.e. 



(k' - 1)/Q' < {k - l/2)/Q < k'/Q'. 



(13) 



This together with integrated Streda formula ( |ll[ ) and 
formulas for the Hall conductance (^,0,^) determines the 
mapping @ and thus the hierarchical tree. For quadratic 
irrationals this algorithm allows to find the tree analyti- 
cally. I do not concentrate here on this matter, but would 
like to make two comments: 

Intermediate fractions. Rational approximants (gener- 
ations of the tree), generated by this procedure, appeared 
to be differ from truncated continues fraction of r\. They 
are: 



V 



O'-i) 



[ni,. . . ,rij - 1] , if a(k) = (-I) 1 1 Qj-i 
[m,...,nj-i], ifcr(k) = (-l) J (Qj -Qj-i) 

(14) 



Thus the parent generation is obtained by either sub- 
tracting 1 from the last quotient n 3 of the continued 
fraction or by truncating the fraction. The sequence of 
generations produced by these iterations is known as in- 
termediate fractions. The golden mean r\ = g -1 is the 
only number which rational approximants are truncated 
continues fraction. 

Takahashi- Suzuki numbers. The entire set of Hall con- 
ductances generated by the flux P/Q are numbers less 
than Q of the form 



These numbers (denominators of intermediate frac- 
tions) are known in integrable models related to U q (SLi2) 
as Takahashi-Suzuki numbers ^0|. They are allowed 
lengths of string solutions of Bethe-Ansatz equations. 
These numbers are also a set of possible dimensions of ir- 
reducible highest weight representations of U q {SLi2) with 
definite parity pi] . 

Each path of the tree may also be characterized by a 



fraction 



,0) 



k/Qj lying on the path and converged 



to a given irrational fraction v. According to eq. (jlJ 
the parent fraction is determined by the daughter one 



as v 



= oh {[Qi-^ - 4j)J + ij with Qi-i = 

Qj — |er(fe)|. The sequence converges to the irrational 
v faster than QJ 1 , i.e., — v\ < cQj 1 . In terms of the 
fractions one may reformulate the scaling hypothesis as 
\Ej — E\ < c\v^ — v\ a ( J \ defining the scaling exponent 
a(J). 

The hierarchical tree, we just described, has been sug- 
gested by the Bethe Ansatz equations for the Harper's 
equation. However, it seems plausible that the construc- 
tion is universal and valid for a general quasiperiodic 
equation, regardless, whether it is integrable or not. A 
set of Hall conductances is the only input of the algo- 
rithm. 

3. Integrability. The Harper's equation (Q) is inte- 
grable as soon as 77 is a rational. Here I adopt a restricted 
definition of the integrability of a linear eqaution: there 
is an isospectral transformation which turns all Bloch so- 
lutions of the Harper's equation to discrete polynomials 
of degree Q. In symbols 



Ak 



m=0 



(16) 



where c nrn is a unitary Q x Q matrix and ^(z) is a poly- 
nomial of the degree Q — 1. In other words, there is a 
gauge (a choice of the gage potential), or a representa- 
tion of the algebra of translations in a magnetic field (2), 
where all wave functions are polynomials. 

In this sence the Harper's equation appears to be in- 
tegrable for any point of the Brillouin zone < k' < 
2n/Q, < k < 2tt JTJ, although the Bethe Ansatz 
equations look especially simple at the so called ratio- 
nal points of the Brillouin zone. The latter correspond 
to the centers and edges of bands. The study of these 
points is sufficient for our purposes. In this case 



Q-i 

j=0 



where aj do not depend on n and p is a constant. 

It appears to be convenient to parameterized polyno- 
mials by its roots: 



3 



Q-i 



0,;Z J 



3=0 



i=0 



(17) 



Below we sketch the results of the Bethe-Ansatz solu- 
tion and skip all aspects of integrability related to cyclic 
representations of U q (Sl2) [p2| . 

Rational points form a zoo. To characterize them we 
introduce parameters r, n, p = ±1. The choice r = 1 
yields levels at the center of bands, while r = — 1 corre- 
sponds to edges of bands. The Chambers relation 



A(fc', k) = detH = 2 cos Qk' + 2A cos Qk 



(18) 



implies that the energy depends on k! and k via A(fc', k). 
Therefore the edges of the energy bands are given by ex- 
trema of A which assumes a minimum/maximum given 
by A = ±(2 + 2A). The middle of bands corresponds to 
A = 0. If P is even (Q is odd) the rational points are la- 
beled by additional discrete parameters k, p — ±1. The 
middle points of bands at k = ±1 are given by the equa- 



tion C ° S j (fc - + ^ 



sin %(k y +n^) 

p 

of bands k' - *k=L=H± 



K-i) 



Q-i 



The edges (r = — 1) 



and k 



l-(-l)~ 
Q \ 2 



21 



Q 2 

are distinguished by parameter p. Being count from 
the bottom of the spectrum, these edges are ordered 
as bottom-top-bottom... if p = — (— 1)t -odd and top- 
bottom-top. . . if p = (— 1)"2 (see [Q) for detailcs). 

For the rational points the transformation ([l6|) is given 
by "quantum dilogarithms" 



m — i 

n( 

3=0 



fc .2n+l/2 A 1/2 + TK P - 1 q-i- 1 / 

X 1 / 2 + Kpqi+V 2 



e q 



(19) 



where p = i exp(i 



). Under this transformation 



the Harper's equation becomes: 

iq(z^ + ztk(— !— (z? - i K ( — )i ) 1 $(qz) 
\ Xqz /V qz 



-iq 1 ( z 2 — itk( 



^)(zi +lK (f)^(zq^) = 
^k\- 1/2 E^(z), (20) 



where one suppose to set z — pq n . However, there is a 
certain advantage to consider the difference equation for 
^(z) in a complex plane z. 

The integrability now reads: all Bloch solutions of the 
difference equation ( |20| ) are polynomials. 

I try to unmasked this transformation by a comment 
bellow, however it becomes meaningful in the U q {sl2) 
setup. 

Let us represent translation operators (^|) by another 
Weyl pair 



UV = qVU, 



(21) 



and setting 



where 



x ~ uv u + b' ± y~ vu u + b' 



a = -ivq 2 A 1/2 
b = -irvq^^X 1 / 2 



(22) 

(23) 
(24) 



Equation (|20|) appears, by choosing a standard represen- 
tation of U and V: 



n+l- 



4. The Bethe Ansatz. Being sure that solutions of the 
Eq.(p0[) are polynomials, we may evaluate it at one of the 
roots of the polynomial Zj. This gives the Bethe-Ansatz 
(BA) equations: 



Q TT 1 z i z k 



q Q E 



k=l 



Zi - qzk 



q? Zi + irnX 2 ) I Zi — inX 



(25) 

Solutions of the BA equations give the wave functions of 
the Harper equation at band's edges and centers. Their 
energy is given by 



E = ipX?q Q {q- q- 1 ) 



A= — tA 



7V2 



7-1/2 



(26) 



The latter is obtained by evaluating the leading terms at 
z — > 00 at eq. j20|) . 

At first glance, the BA equations ( p5| ) look even more 
complicated than the original Harper equation. This is 
true as long as Q is not large. However, the BA equations 
( p5| ) provide a better description of the problem in the 
most interesting, incommensurate, limit P,Q~^ 00, r\ — > 
irrational number. 

At P odd the BA equations admit an exact zero mode 
solution at for E = p3] |. For a quasiclassical analysis 
of the BA equations at 77 — > 0, see p4| . 

Below, we consider the most interesting case A = 1. 

3. String hypothesis. Here we formulate the string 
hypothesis which allows us to obtain the solutions of the 
BA equations (at A = 1) with an accuracy 0{Q~ 2 ). The 
hypothesis is based on the analysis of singularities of the 
BA, and is supported by extensive numerics jj. Here 
we just formulate the string hypothesis and present some 
immediate consequences. To proceed, we need the notion 
of strings. 

A string of spin I , parity v— ± 1 and center xi is a set 
of 21 + 1 complex numbers: 



-L-l + l. 



J. 



(27) 



which have a common modulus x\ > (a center of the 
string), a parity vi = ±1 and differ by multiples of qi. 

Now we are ready to formulate the string hypothesis 
— a central concept of this analysis: 
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• At large Q each solution of the BA consists of 
strings. 

• Each solution can be labeled by spins {lj, ■ ■ ■} 
and parities {vj,Vj-i, . . .} of strings, such that the 
total number of roots + 1) = <9 — 1- We 
refer to the set of lengths and parities of strings 
constituting the solution for a given energy level as 
to a string content of this level. Not more than 
two strings with a given length and parity may be 
found in a string content of the solution. 

• The length of the longest string in a string content 
of a given energy level is the Hall conductance of the 
corresponding band: 21 + 1 = |u(fc)|. The period of 
this string qi — e' l7 ' ri1 is uniquely determined by the 
requirement that r\i = is the best approximant 



for the period r/, so that q\ 



2;+i 



±1. 



• The parity of the longest string is vi = —iql +1 ^ 2 i> — 
(—1)^1/. The center of the longest string is xi — 

1 + 0(1/0. 

• The remaining roots of the state is a solution of 
the BA equation for the parent state of the parent 
generation. 

The string hypothesis states that 
l 

^daugther^j ^ JJ ( z _ a^p)** 8 ™"*^). (28) 
m=— I 

and that 21 + 1 = cr(fc)| is the absolute value of the Hall 
conductance of the "daughter" band. 

This simple hypothesis allows one to construct a com- 
plete set of wave functions by virtue of the iterative pro- 
cedure: 

Starting from an irrational 77, we first generate a hier- 
archical tree. Let us choose a branch of the tree J. Find 
the Hall conductance of the bands belong to the branch 
down to the origin. This determines the lengths, periods 
and parities of the string content and therefore zeros of 
the wave function of the chosen branch. 

The only unknowns are the centers of strings. They, 
however approach 1 with accuracy 0(1/1). The accuracy 
of the recursive eq.(28) is 0(l~ 2 ). The string content of a 
state (i.e., lengths and parities of strings) is a topological 
characteristics, while the centers of strings are not. 

The strings hierarchy has been obtained by analysis of 
singularities of BA (see || ) . As it was expected a set of 
possible lengths of strings is a set of Takahashi-Suzuki 
numbers, known in the Bcthc- Ansatz literature. Eq.(0) 
provides a relation between them and Hall conductances. 

To illustrate the iterative procedure let us consider the 
bottom edges of the lowest band of the spectrum and 
choose r\ = to be the golden mean. The sequence 

of rational approximants is given by ratios of subsequent 



Fibonacci numbers rji — — pr 1 , where the Fi are Fibonacci 
numbers (Fi = F;_ 2 + iCi and F = F x = 1). The 
set of Hall conductances = Takahashi-Suzuki numbers = 
allowed lengths of strings are again Fibonacci numbers: 
Qk-i = ffc-i- The considered branch of hierarchical tree 
connects edges (r = — l,/t=l,/i=l) of the lowest bands 
of generations r/^k = i^fe-i/i^fe. Their string content 
consists of pairs of strings with length 2l n + 1 = F^ n+ i, 
n = 0, 1, . . . , k — 1, parities +1 and inverted centers Xk 
and xZ . According to the string hypothesis the wave 
function of this state is 

fc-l In 

*(*i?73fc) « n n ( z - x ijji z - x i!4j- ( 29 ) 

n=0 j=-l n 

Centers of the strings xi n are close to 1 but can not be 
obtained from the string hypothesis alone. 

5. Gaps. A direct application of the string hypothesis, 
suggested by J. Bellissard, is the calculation of the gap 
distribution p(D) 1 i.e., the number of gaps with magni- 
tude between D and D+dD. The result is p(T>) ~ X>~ 3 / 2 , 
which essentially means that the width of the smallest 
gaps scales as D m i n ~ 1/Q 2 . This result confirms nu- 
merical analysis of Ref . (25| . 

6. Scaling hypothesis and finite size corrections: stat- 
ing the problem. The string hypothesis solves the Bethe 
Ansatz equations with an accuracy 0(Q~ 2 ), i.e., is 
asymptotically exact in the incommensurate limit Q — > 
00. It alone gives the explicit asymptotically exact form 
of wave functions and provides the hierarchical tree and 
topology of the Cantor set spectrum. However, the most 
interesting quantitative characteristics of the spectrum 
are hidden in the finite size corrections of the order of 
Q~ 2 to the bare value of strings. Among them are the 
anomalous dimensions of the spectrum e j. They depend 
on the branch and on arithmetics of rj (according to ref. 
|p6| exponents ej vary between 0.171 and —0.374 for the 
golden mean 77 = v/ ^~ 1 ). Can anomalous dimensions be 
found analytically, by finding finite size corrections to the 
string solutions? This is a technically involved but a fas- 
cinating and important problem. Its solution would pro- 
vide an ultimate information of the spectrum and most 
interesting physical properties of the system. It also may 
suggest the conformal bootstrap and operator algebra ap- 
proach, which has been proven to be effective for finding 
the finite size corrections of integrable systems, without 
the actual solving the Bethe Ansatz. 

7. I would like to acknowledge inspiring numerics pro- 
vided by Y. Hatsugai during the initial stages of formu- 
lated the string hypothesis, and useful discussions with 
Y. Avron, G. Huber, M. Kohmoto, R. Seiler and A. 
Zabrodin. I also would like to thank M. Ninomiya for 
the hospitality I received in YITP, where the text of this 
lecture was written. 
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